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Abstract 

Previous work has shown that Benjamin-Feir unstable traveling waves of the com- 
plex Ginzburg-Landau equation (CGLE) in two spatial dimensions cannot be sta- 
bilized using a particular time-delayed feedback control mechanism known as 'time- 
delay autosynchronisation'. In this paper, we show that the addition of similar 
spatial feedback terms can be used to stabilize such waves. This type of feedback 
is a generalization of the time-delay method of Pyragas (Phys. Letts. A 170, 1992) 
and has been previously used to stabilize waves in the one-dimensional CGLE by 
Montgomery and Silber {Nonlinearity 17, 2004). We consider two cases in which 
the feedback contains either one or two spatial terms. We focus on how the spatial 
terms may be chosen to select the direction of travel of the plane waves. Numerical 
linear stability calculations demonstrate the results of our analysis. 



1 Introduction 



During the past two decades considerable progress has been made in our un- 
derstanding of the spontaneous emergence of patterns in spatially-extended 
non-equilibrium systems. The existence of simple spatial or spatio-temporal 
patterns has been rigorously established on the basis of equivariant bifurca- 
tion theory (see [T] and references therein). However, these simple patterns are 
often unstable in a given system, which evolves instead to a state of spatio- 
temporal chaos. Patterns that result from a symmetry-breaking Hopf bifur- 
cation seem to be especially vulnerable to instability. A current challenge in 
pattern-formation research is to develop control schemes that stabilize the 
simple patterned states, so that a desired, otherwise unstable, solution may 
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be realized. This paper extends current research on feedback control of oscil- 
latory patterns, focusing on traveling wave solutions of the two-dimensional 
complex Ginzburg-Landau equation. One of the goals of this line of research 
is to develop control algorithms that exploit the underlying symmetries of a 
targeted pattern in order to stabilize it in a non-invasive fashion. That is, 
if the pattern is invariant under some spatial, temporal or spatio-temporal 
transformation, then it may be possible to construct a feedback based on this 
transformation which vanishes when the target pattern is achieved. For exam- 
ple, if the pattern has the property that A(x, t) = A(Rx, t — r) where i?x is 
some Euclidean transformation of the spatial variable, and r > is a possible 
time delay, then a feedback which is proportional to A(i?x, t — r) — A(pc, t) 
will be non-invasive. 

Ott, Grebogi and Yorke [2] were the first to develop a control algorithm that 
suppresses chaos in favor of a simple unstable periodic orbit (UPO) in low- 
degree-of-freedom chaotic systems. In this approach, small perturbations are 
applied to a system parameter in order to keep the system close to the UPO. 
However, this method requires constant monitoring of the system, and also 
may not be effective in rapidly evolving systems. 

Pyragas [3] introduced a second approach, sometimes called 'time-delayed au- 
tosynchronisation' (TDAS), in which the feedback is proportional to the differ- 
ence between the current and a past state of the system. That is, the feedback 
is F = j(A(t) — A(t — At)) where At is the period of the desired UPO, and 
A(t) is some state variable. This method has a number of attractive prop- 
erties. First, the UPO of the original system will also be a solution of the 
system with feedback, and so control may be achieved in a non-invasive man- 
ner. Second, the only information required a priori is the period of the desired 
UPO. The method has been implemented successfully in a variety of labora- 
tory experiments on electronic [HE], laser [6], plasma [7f8f9] . mechanical [ID] 
and chemical systems [TT|12jll3lfT3] ; more examples can be found in a recent 
review by Pyragas [15]. 

An extension of TDAS proposed by Socolar et al. [H] incorporates informa- 
tion at many previous times and is known as 'extended time-delayed autosyn- 
chronisation' (ETDAS). This method was used by Bleich and Socolar [17] to 
stabilize unstable traveling wave solutions of the complex Ginzburg-Landau 
equation (CGLE) in one spatial dimension. For spatially-extended pattern- 
forming systems, other proposed modifications of the time-delay autosynchro- 
nization scheme of Pyragas include global feedback control [T81I19II20] . where 
the magnitude of the feedback depends on some spatial average, or maximum, 
of a quantity at an earlier time. Yet others take into account the spatial peri- 
odicity, as well as temporal periodicity, of the targeted pattern. For instance, 
Lu, Yu and Harrison [22] used numerical simulations of the two-dimensional 
Maxwell-Bloch equations describing a three-level laser system to demonstrate 
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that spatio-temporal chaos could be suppressed by applying a linear combi- 
nation of time-delay feedback and an analogous spatially-translated feedback 
term of the form: 

F s = p{[E(x + x , y, t) - E(x, y, £)] + [E(x, y + y , t) - E(x, y, t)}}. 

Here di = (x , 0) and d 2 = (0, y ) are translation vectors associated with the 
feedback; the feedback is non-invasive when the resulting pattern is periodic 
in each of the dj-directions, with spatial periods \dj\. Their numerical investi- 
gations reveal that these two spatial controls somehow select the direction of 
travel of the stabilized plane wave. This type of feedback, with a combination 
of temporal and spatial terms, was extensively studied by Montgomery and 
Silber [21] in the context of the one- dimensional CGLE. 

The aim of our paper is to extend the results of [21J to the CGLE in two spatial 
dimensions. This problem is interesting since Harrington and Socolar [23J have 
shown that unstable traveling waves in the 2D CGLE cannot be stabilized 
using only temporal feedback, due to the presence of torsion-free modes (that 
is, modes which have purely real Floquet multipliers). We circumvent this 
difficulty by using a combination of temporal and spatial terms, as in |21j, and 
additionally consider the effect of having either one or two spatially shifted 
terms in order to steer the direction of the traveling waves in the plane. 

Our analysis is based on a linear stability analysis of the traveling wave solu- 
tions of the CGLE. The analysis leads to a system of delay differential equa- 
tions (see [2¥|2"5] for more on DDEs). We find stability boundaries by searching 
for critical curves/surfaces of Hopf bifurcations within the system of DDEs. 
The analysis is made possible because the CGLE admits an exact family of 
solutions in the form of a traveling plane waves RQ^^+^t ) parameterized by 
the wave vector k. For this problem the dispersion relation that relates |k| 
and uj is known precisely. 

We expect that our methods could equally well be applied to other spatially 
extended systems for which traveling plane waves exist. Such solutions arise 
for instance when a spatially-uniform solution undergoes a Hopf bifurcation. 
However, in many instances, the form of the wave and the dispersion relation 
are unknown, or known only approximately. In this case, we propose that 
just one spatial term be incorporated in the feedback, and consider target 
waves that travel in a direction relative to this feedback term so that both the 
temporal and spatial feedback terms vanish. In this way the wavenumber and 
frequency of the plane wave can be matched even when the dispersion relation 
is unknown, although the direction of travel must be left free. 

This paper is organized as follows. In Section [2] we review some properties of 
the CGLE and the feedback we are applying, and in Section [3] we describe 
previously known results regarding these systems. In Section @] we discuss 
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the effects of additional spatial feedback and how the spatial shifts must be 
chosen in order to stabilize the wave. Section [5] contains a geometric way to 
think about the results of [2TJ and in Section Owe provide a numerical example 
to demonstrate our results. Section [7] concludes. 



2 Problem setup 

2.1 The complex Ginzburg-Landau equation 



In this section we review some properties of traveling wave solutions of the 
complex Ginzburg-Landau equation (CGLE). The CGLE is an amplitude 
equation for some disturbance A(x,t) in a spatially extended system near the 
onset of a Hopf bifurcation that sets in with zero wavenumber, i.e. a spatially- 
uniform oscillatory instability. After appropriate rescalings and in a frame that 
is oscillating at the Hopf frequency, it can be written generically as: 

r)A 

_ = A + (1 + i&O V 2 A - (b 3 - i)\A\ 2 A, (1) 

where bi and b 3 are real parameters. We consider only the situation just after 
a supercritical Hopf bifurcation, so the coefficient of the linear term is positive 
(and has been rescaled to equal unity) and the parameter b 3 is assumed to be 
positive. In this paper, we consider the spatial extent of the problem to be 
two-dimensional, so x = (x,y). We consider traveling wave solutions of the 
form 

A TW = Re lk - M , (2) 

where the amplitude R and frequency u> are determined by the wavenumber 
k = |k| through the relations 

i? 2 = l^!, uo = R 2 -b 1 k 2 . (3) 

Since 63 > it follows from ([3]) that k < 1. We also assume that k 7^ 0. 

In the Benjamin-Feir unstable regime (b\ > 63 > 0) all solutions of the form (j2j) 
are unstable. In this paper we investigate how the addition of feedback terms 
to (pQ) affects the linear stability of these solutions in this regime. Specifically 
we consider 

f) A 

_ = A + (1 + ih)V 2 A - (63 - i)\A\ 2 A + F, (4) 
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where F is a feedback term given by 

N 

F = j[A{x,t) -A(x,t- At)} +^ Pj [A(xL + Axj,t) - A(x,t)]. (5) 

i=i 

The delay parameter At is positive, and we restrict our analysis to the case 
where 7 (the gain) and pj are real. We consider the two cases N = 1 and 
N = 2, that is, the feedback F can contain either one or two spatially-shifted 
feedback terms. 

In this paper, we study the stability of traveling waves for which the feedback 
terms vanish. We choose the spatial and temporal shifts Ax,,- and At such that 

At\u\ = 2tt Ax,- ■ k = 27T7ij, rij eZ, j = 1, . . . , N (6) 

for a specific targeted traveling wave of the form (121), with wavevector k = 
(fca., and frequency to. 

There may be many other solutions to (jl]) for which the feedback term does 
not vanish. We do not study these solutions in this paper - our analysis is a 
purely local stability analysis assuming we are already close to the targeted 
wave. 

Since we know the dispersion relation (|3J) relating the wave frequency u and 
the wavenumber k = |k|, the spatial and temporal shifts are related and cannot 
be chosen independently. We now consider how the shifts must be chosen in 
the two cases N = 2 and N = 1. 

Note that if we were to choose N > 2, there would be some relations between 
the Axj in order that they satisfy (jH]). For instance, for N = 3, we would need 
Ax 3 = rii Ax 2 + n 2 Ax 2 for some integers n\ and ri2 ■ This is because we are in 
two spatial dimensions. We do not consider the cases of higher N here, but it 
may be that using additional shifted terms can increase the region of stability 
of the traveling waves, much in the same way as ETDAS [17J. 

For iV = 2, we insist that Axi and Ax 2 are not parallel (Axi x Ax 2 7^ 0). 
Then for specified Ax x and Ax 2 there is a two-dimensional dual lattice of 
possible target wavevectors k. The lattice generators kx and k 2 satisfy 

^ = lAx'xlLr VA* 3 = 0, j = 1,2. (7) 

The frequency of the desired wave is specified by the choice of At, which 
must be consistent with the choice of k, since the frequency and wavenumber 
are related a priori by the dispersion relation ([3]). (This 'overspecification' of 
the targeted wave is analogous to the one-dimensional case of [2T], where one 
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temporal and one spatial term select u and k which are already related by the 
dispersion relation.) Note that if k lies on the lattice, so also does — k, and 
thus the direction (right/left) of travel of the wave is not selected. 

If the same procedure were repeated for a different amplitude equation for 
which the dispersion relation were not known, then in the 'overspecified' N = 2 
case there could be a discrepancy between the temporal and spatial feedback 
terms, and it may not be possible for both to vanish. However, in the case 
N — 1, there is more freedom, as we explain below. 

For AT = 1 it is not possible to choose Axi such that the spatially-shifted term 
targets a single traveling wave. That is, for a given ni, there is a continuum 
of wavevectors k for which k • Axi = 27rni and the spatial feedback term 
vanishes. The possible target k for a given Axj are indicated by the bold lines 
in Figured], in section HI However, the choice of At selects some u, and hence 
|k| is also selected, by the dispersion relation. The direction of k (again, up to 
a sign) is selected by the spatially shifted term, and so a single traveling wave 
can be targeted. This is explained in more detail in Section 14.2.11 

Note that F can only improve the stability of the traveling waves if 7 < and 
likewise pj > 0. The reasons for this are explained in detail in [2T] and we will 
only consider this parameter regime. 



2.2 Stability analysis 



The linear stability of traveling waves (T5]) is calculated by considering the 
effect of small amplitude perturbations, for some perturbation wavevector q = 
{q x ,Q y )- We write 



A = Re 1 * 



l + a+(t) e iq - x +a_(t) e" 



iq-x\ 



substitute into (jl]) and linearize in a + and a_. This results in a system of 
delay differential equations for a + and the complex conjugate of a_ (denoted 
by a!.): 






a+(t - At) 
al(t-At) 



where 



J 



- Cl q 2 - c 2 R 2 -c 2 R 2 

-c\q 2 - c\R 2 



-c* 2 R 2 



+2kq 




iq- Ax 



(9) 



1 
(10) 
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with c\ = 1 + ibi, C2 = 



In the following sections we review previous results concerning feedback control 
in equations similar to (jUJ) and then present our new results. 



3 Previous results 



Our stability analysis builds on results of Harrington and Socolar [23], Naka- 
jima [26J and Montgomery and Silber [2T]. For completeness, in this section we 
give a summary of the results of [23] and [26] , and also of previously known re- 
sults concerning some instabilities of the CGLE without feedback. We discuss 
details of [21] in Section [5j 



3. 1 Instabilities with no feedback terms 




In the absence of any feedback terms (i.e. 7 = pj = 0), traveling waves are sta- 
ble to perturbations with wavevector q if both eigenvalues of the Jacobian (fit)]) 
have negative real part. 

The Benjamin-Feir (B-F) instability is a long wavelength instability, that is, 
perturbations with q — |q| sufficiently small will grow. In the limit q ^ 1, 
a simple calculation shows that one eigenvalue of J always has negative real 
part, and the other has real part equal to: 

(k-q) 2 + 0(g 3 ). 

Hence, if bi > b 3 > (the B-F unstable regime), then this is positive for all 
k, and long wavelength perturbations will grow. 

A second type of instability important to our subsequent analysis, is associated 
with perturbations with k • q = (q no longer necessarily small). It can be 
shown [23] that these perturbations grow for small enough q, that is, those 
with 

q < & = x + 62 ( n ) 
Perturbations with k • q = and q > q cr decay. 

When q is sufficiently large, the eigenvalues of J ~ — q 2 , so all sufficiently 
shortwave perturbations will decay. However, there may be other regions of 
unstable q in addition to those described above. 
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3.2 A result of Harrington and Socolar 

We now summarize a result from Harrington and Socolar [23] concerning the 
stabilization of plane waves in the 2D CGLE with only temporal feedback. 
Their result relies on a proof of Nakajima |26j which has recently [2"T1I2"%] been 
shown to be incorrect. However, as we discuss below, the result of Nakajima 
does still hold in some cases, and the results of [23] are still correct. 

The claim of Nakajima is as follows. Suppose an ODE 

x = f(x{t),t), xeW 1 , (12) 

contains an unstable periodic orbit, x*(t), with period T. Feedback is added 
to (I12p so it is now of the form 

x (t) = f{x{t),t) + K{x{t) - x(t - T)), (13) 

where K is an n x n gain matrix. Then if the linearization of (|T2|) about the 
orbit x*(t) has an odd number of positive, real unstable Floquet multipliers 
then there is no value of K for which x*(t) is a stable solution of (1131) . 

The flaw in the proof of Nakajima is that the neutral Floquet multiplier asso- 
ciated with perturbations in the direction of the periodic orbit was neglected. 

However, in both the example of Harrington and Socolar, and in our work, 
the analysis of the stability of the traveling wave has been reduced to the 
study of the stability of a fixed point (namely, the zero solution to the linear 
DDE ([2])). The result of Nakajima can be applied to hyperbolic fixed points, 
since they do not have a trivial Floquet multiplier, and so the results of [23] still 
stand. (A similar result for fixed points in maps, rather than flows, is given 
by Ushio [29].) In our examples, the perturbation along the periodic orbit 
which yields a neutral Floquet multiplier is associated with a perturbation 
with wavevector q = 0. It is not this type of perturbation we are concerned 
with in the following. 

The result of Harrington and Socolar is summarized as follows. Consider Q 
with no spatial feedback, that is, pj = 0. In 2D there always exists an unstable 
perturbation wavenumber q which satisfies k ■ q = 0. In this situation, the 
Jacobian J is equal to 

'-c x q 2 -c 2 R 2 -c 2 R 2 \ 
K -c* 2 R 2 -c\q 2 - c\R 2 ) 

where Re(ci), Re(c2) > 0. If < q < q cr , then J has real trace and real, 
negative determinant, hence has real eigenvalues of opposite sign. 
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Note that in only one spatial dimension, we can only have k • q = if either 
q = or k = 0, so the same result does not apply. The 'flat' perturbations 
with q = do not cause instabilities for 7 < 0, which is the parameter range 
we are interested in. 



4 Additional spatial feedback 

In this section we extend the results described above to include the effects 
of additional spatial feedback, with either one or two spatial terms. This sec- 
tion is organized as follows. First, we show that there are some perturbation 
wavevectors q which are unaffected by the spatial feedback. We then discuss 
how the displacements Axj must be chosen to target a particular wavevector 
k, for both the cases N = 1 and N = 2. 

4-1 Perturbations unaffected by spatial feedback 

Consider the Jacobian matrix J given in (fTUj) . If the perturbation wavevector 
q satisfies q • Axj = 2irmj, for some rrij G Z, then J is the same as if pj = 0. 
Hence, these perturbations q are not affected by the spatial feedback. Com- 
bining this with the results of Harrington and Socolar [23] described in the 
previous section, if these q are unstable for the system with no feedback, and 
additionally satisfy k - q = 0, then the traveling waves cannot be stabilized for 
any values of pj and 7. That is, if there is a q 7^ which satisfies 

k ■ q = 0, q ■ Axj = 2nnj, rij G Z and |q| < q cr , (14) 

then the traveling waves cannot be stabilized. 

If there are no such q, then it is possible that there exists some choice of pj 
and 7 for which the traveling wave is stable. 

4-2 Choice o/Axj 

Recall that in order for the spatial feedback to vanish at the targeted wave 
solution, we require that 

k ■ Axj = 27mj, rij G Z. (15) 

In one spatial dimension, this gives a unique choice of Axj (for say, rij = 1) for 
a given k, but in two spatial dimensions this is not the case. We now discuss 
this further. 
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Fig. 1. The bold lines indicate those points in k-space which satisfy k • Axj = 2-7rni 
(m € Z). The vector k p is parallel to Axi. The circle indicates the wavevectors k for 
which the temporal feedback term vanishes with |u;|At = 2tt. The intersections of 
the circle and the lines (marked by dots) are the possible target k. As the direction 
of k gets closer to k p , either by changing the size of the circle, or altering the angle 
of Axi, the magnitude of the q which satisfy k • q = and also lie on the bold lines 
increases. 

4.2.1 The case N=l 

For N — 1 there is a continuum of wavevectors k for which k • Axi = 2nri\ 
and the feedback term vanishes. The possible k are indicated by the bold 
straight lines in Figure [TJ However, isolated k can be picked out from this 
continuum by the temporal feedback term (which fixes At and hence u) and 
the dispersion relation. The dispersion relation fixes k 2 for a given u and so 
gives a circle of possible wavevectors. If we let |o;|At = Inn for some integer 
n > 1, then this will give a family of circles, of decreasing k as n increases 
(since k must decrease as uj increases), however, we only consider n = 1. The 
intersection of this circle and the straight lines in Figure Q] give the possible k 
for which the feedback term vanishes. In this example, there are three possible 
k (up to a sign). As we explain below, the k with direction closest to that of 
Axi will be the easiest to stabilize. However, it is certainly possible that more 
than one of these k may be simultaneously stabilized. 

If the dispersion relation for a system is known, as is the case for the CGLE, 
then the location of this circle is known, and At and Axi can be chosen to 
target a particular traveling wave. If the dispersion relation is not known, 
then the location of the circle is not known, but it still exists. In this case, 
we can, say, choose At to pick the frequency of the target wave, and also pick 
some Axi. We will not be able to choose in advance the wavevector k of the 
resulting targeted wave, but the possible k for which the feedback vanishes will 
be isolated in k-space. This method provides us with a possible mechanism for 
using this stabilization scheme for a system in which the dispersion relation is 
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not known. 



If there is a q which satisfies (I14p . it will also lie on the bold lines in Figure [TJ 
However, for a given k, Axi can be chosen in such a way that there are no q 
which satisfy the resonance conditions (1141) . To see this, consider the situation 
where Axi is chosen to be almost parallel to k, then the only q perpendicular 
to k which also lie on the bold lines will have a very large magnitude and 
hence do not correspond to wavenumbers of destabilizing perturbations. Once 
Axi is chosen to avoid the resonance conditions (j!4p the methods developed 
in [2T] can then be used to determine the stability of the traveling waves. 

By analogy with the one-dimensional case, we might hope that we are able 
to stabilize the wave for which k and Axi are parallel (indicated by k p on 
Figure [T]). However, this is not the case since there always exist unstable 
longwave perturbations (q small) which satisfy k p • q = q ■ Axi = 0, and these 
perturbations are unaffected by the spatial feedback. Hence this plane wave 
cannot be stabilized in this fashion. 

4.2.2 The case N = 2 

For the spatial feedback with two terms to vanish at the targeted wave solu- 
tion, k must satisfy 

k ■ Axi = 27rni, k • Ax2 = 2irn2, n\, rii G Z. 

This gives a lattice of possible target wavevectors with generators as given 
in equation ([7j). An example lattice is shown in Figure [21 Recall that for the 
CGLE we must have k — |k| < 1, so for a given target k, it is possible to choose 
Axi and Ax2 such that k and — k are the only non-trivial lattice points inside 
the unit circle. However, for a specific choice of k, the choice of Axi and Ax2 
resulting in an appropriate lattice is not unique. 

The perturbations q which satisfy (|14[) will lie on the same lattice as k. If q cr 
is small enough, there will be no unstable q (that is, with |q| < q cr ) which lie 
on the lattice. Therefore, so long as Axi and Ax2 are chosen appropriately, 
the theorem of Nakajima does not apply. The methods developed in [21] can 
then be used to determine the regions of stability for these waves. In the next 
section we review the approach taken in [21]. 



5 Analysis 

In this section we discuss a result found in Montgomery and Silber [21]. They 
show that for the CGLE in ID with feedback as in (j3J), the 'best' choice of 7 is 
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Fig. 2. The figure shows a lattice of possible k-vectors for which the spatial feedback 
vanishes in the case N = 2. The bold lines indicate the lattice generators. The 
shaded region shows those k with |k| < 1. In this example only two non-trivial 
lattice points lie within this disc. 

— 1/ At. That is, if it is possible to stabilize a traveling wave for some choice of 
feedback parameters, this can certainly be achieved by choosing 7 = — I /At. 

The proof given in [21] can be extended to two dimensions. Here we omit 
those details and give a briefer, and more intuitive reason as to why their 
result holds, although it is not rigorous. 

We first reduce the delay differential equation Q to a single equation by 
diagonalizing the matrix J to give 

a(t) = mia(t) +j{a{t) - a(t - At)) (16) 
b(t) = m 2 b(t)+7(b(t)-b(t-At)) (17) 

where a and b are appropriate linear combinations of a+ and a*_, and 

2 

m k = m k + ]T Pj (e^ A ^ -1), k = 1, 2. 
i=i 

The rrtk are the eigenvalues of the Jacobian ([TO]) with no feedback (that is, 
Pj = 7 = 0), with Re(mi) > Re(m 2 ). As described in [21J, the b equation 
does not affect the stability boundary since it is associated with the more 
stable eigenvalue. We therefore concentrate on the a equation, and write ih\ = 
a(q) +z/3(q), where a and (5 are real functions of q. Hence, a(q) is the largest 
real part of the eigenvalues of J when there is no temporal feedback. Since the 
feedback can only improve the stability of the traveling waves, instabilities are 
only possible when a(q) > 0. 
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Fig. 3. Schematic diagram showing critical curves in unstable (left-hand picture) 
and stable (right-hand picture) cases. The transition between the two takes place 
at 7 = — 1/At. The regions of the diagram are marked as containing perturbations 
q which are stable (s) or unstable (u). 

Critical surfaces describing the boundaries between regions of stable and un- 
stable perturbations have a(t) = e xt with A = iu, v G R (i.e. they are points 
of Hopf bifurcations since v ^ 0). We can write equations for these surfaces 
in 7-q space, parameterized by v, as: 

ct(q) = 7(cos v — 1) (18) 
7 At sin z/ = v + At/?(q). (19) 

The results of [2T] can be used to give an algebraic condition for the stability 
of the traveling waves. That is, if there are no solutions u, q to the equations 

cos v - 1 = -Ator(q) (20) 
sin v + v = -At/3 (q) (21) 

then the traveling wave can be stabilized at 7 = —1/At. 

In Figure [3] we give a sketch of the shape of the curves given by ({TBI and (|T9l) 
for a fixed direction of q, in two cases, one before and one after a change of 
stability of the traveling waves. At points with cv(q) = 0, the curves either 
intersect 7 = (if v 7^ 27m), or asymptote to 7 = ±00 (when v = 27m). 

These curves divide regions of stable and unstable perturbations. For a travel- 
ing wave to be stable, it must be stable for all perturbations q. One mechanism 
by which the stability of a traveling wave can change is for two critical curves 
to collide at points with |^ = 0. Figure [3] shows a schematic of this transition. 
Of course, the stability of the wave could also hypothetically change without 
the two curves in the right hand picture of Figure [3] colliding; the minimum 
of the top curve would just have to extend below the maximum of the lower 
curve. We do not consider this mechanism of stability change here, but it is 
covered by the analysis in [21] . 
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Hence, we wish to identify points on the critical curves which have = 0, 
since it is at these points that the surfaces can collide and change the stability 
of the traveling waves. Differentiating (ITS]) and ffT9]) with respect to 7, and 
setting |s = gives: 

1 (22) 
7 At cos z/) (23) 

Eliminating results in 

(cos v — 1)(1 + 7At) = 0. 

When cosz/ = 1, the critical curves asymptote to 7 = ±00, a(q) = 0, so 
the only turning points with ^ = will be at 7 = — which is the value 
identified as optimal in [21] . 



av 

7 sin i/—— = cos v — 
c>7 

At sin z/ = — ( 1 — 

c/7 



6 Numerical results 



We now give two numerical examples to show that there are regions of pa- 
rameter space in which traveling waves can be stabilized using the methods 
described above. The first example we give can be stabilized for N = 2 but not, 
for the parameters we use, for N = 1. The second example can be stabilized 
for N = 1. We also demonstrate the way in which stability of the traveling 
wave is lost, as described above - two critical surfaces collide at 7 = —1/At. 

We used the Matlab package DDE-BIFTOOL [30] to analyze the linearized 
system (Q and search for Hopf bifurcations in order to locate the critical 
surfaces between stable and unstable perturbations in 7-q space. For our first 
example, the parameters used are b\ = 2.5, 63 = 2. Without loss of generality, 
we choose k to be parallel to the x-axis. The spatial shifts Axi and Ax2 are 
chosen to be parallel to the vector (1,1) and k respectively, and also to satisfy 
Ax,- ■ k = 2tt, j = 1, 2. We choose pi = 0.01 and p 2 = 0.007. Recall that the 
wave cannot be stabilized if there exist unstable perturbations with k • q = 
and q • Axj = 27mj. In this example, q cr = 2( - 1 29 fc ^ and it is simple to check in 
each case that there are no such q with |q| < q cr . 

The contour plots of a(q) in Figure H] show the regions of q-space in which it is 
possible to find instabilities - that is, those regions where a(q) > 0. Note that 
there is a symmetry q — > — q. In the case N = 1 there are two disjoint regions, 
but for N = 2 there is only one. If there were no temporal feedback (i.e. 7 = 0), 
then the traveling wave would be unstable. In this example it is not possible 
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(c) N = 2 (d) N = 2, pi = 0.1, p 2 = 0.07 

Fig. 4. The figures show contour plots of a(q) close to q = in the three cases 
N = 0, 1,2. As |q| — > oo, a(q) — > — oo. Note that a(q) corresponds to the largest 
real part of the eigenvalues of the linearized system with no temporal feedback (i.e. 
7 = 0). The regions which have positive eigenvalues (corresponding to instabilities) 
are shaded in. Parameters used are b\ = 2.5, 63 = 2, and k = (0.285,0). Axi is 
parallel to (1, 1), Ax2 is parallel to k, and Ax,- k = 2tt. In (a), (b) and (c), p\ = 0.01 
and P2 = 0.007, which are the parameters used in later investigations. In (d) p\ and 
P2 are increased by a factor of ten, and there are still unstable regions, which cannot 
be made to vanish no matter how large we choose the pj. The plots of a(q) for the 
other values of k used in Section [6] are very similar, and so we do not show them 
here. 



to suppress all the instabilities using just the spatial feedback; figure IU(d) has 
spatial feedback ten times that of (a)-(c). The q which are still unstable are 
centered around points which satisfy q • Axj = 2nm for some m G N, so 
the terms proportional to pj in ( flOl) vanish (i.e. these q are unaffected by the 
spatial feedback). Therefore, however large the pi are in this example, these q 
will still be unstable. It may, in other examples, be possible that the traveling 
wave can be stabilized using only spatial feedback. Examples of this for the 
one-dimensional CGLE are shown in [21]. We now give more details on the 
addition of temporal feedback to the two cases N — 1 and N = 2. 
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Qx Hv q y qx 

(a) |k| = 0.285 (b) |k| = 0.287 

Fig. 5. The two plots show surfaces of Hopf bifurcations for two values of |k|. The 
regions above the upper surface and below the lower surfaces are regions of unstable 
perturbations. In (a) the surfaces are disjoint and so there is some range of 7 for 
which all perturbations will decay and so the traveling waves can be stabilized. In 
(b) this is not the case. Notice that the range of q in these plots includes that where 
a(q) > (see figure ID^c)). Compare also with figure [H which shows plots of the 
Floquet multipliers at 7 = — 1/Ai. Parameters used are b\ = 2.5, 63 = 2, p\ = 0.01, 
P2 = 0.007. Axi is parallel to (1,0) and Ax2 is parallel to (1, 1). 

6.1 N=2 

Parameters are chosen as described above. We consider two values of k, k = 
(0.285, 0) and k = (0.287, 0). In both cases, the value of -1/At « -0.04. 

Figure [5] shows bifurcation surfaces for the two cases. In the first case, the 
surfaces are disjoint and so the traveling wave can be stabilized for a range of 
choices of 7, which include —1/At. In the second case, the two surfaces have 
joined, and there is no choice of 7 for which all q perturbations are stable. 

Figure [6] shows the real part of the Floquet exponents of the linearized system 
at 7 = —1/At. We can see that for |k| = 0.285, the growth rates are nega- 
tive for all q, so the traveling wave is stable. For |k| = 0.287, there are two 
(symmetry- related) regions in the q-plane which have positive growth rates 
and hence the traveling wave is unstable. 

6.2 N=l 

We now investigate whether the traveling waves can be stabilized using only 
one spatially shifted term. We first investigate the same example as used above 
for N = 2, but only use the spatial shift Axi which is parallel to the vector 
(1,1). (Note that we cannot use only the second spatial shift Ax 2 since it is 
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(c) |k| = 0.283, N = 1 (d) |k| = 0.285, N = 1 



Fig. 6. The contour plots show the real parts of the Floquet exponents of the lin- 
earized system at 7 = —1/At ~ —0.04. The upper two plots show |k| = 0.285, 0.287, 
N = 2, and the lower two show |k| = 0.283, 0.285, N = 1. The regions with positive 
growth rate are shaded in, as in Figure [H We can see that for N = 2, the case 
|k| = 0.285 is stable but the case |k| = 0.287 is unstable. For N = 1 there are 
additional regions of instability so the wave is not stabilized. Parameters used are 
b\ = 2.5, 63 = 2, pi = 0.01, P2 = 0.007. Axi is parallel to (1, 0) and Ax2 is parallel 
to (1,1). 

parallel to k.) Figure [6] shows the real part of the Floquet exponents of the 
linearized system at 7 = —1/At for this case, for two values of |k|. Note that 
there are more regions of instability in the N = 1 case than the N = 2 case. 
The instability at q ~ (0.3, 0) is present in both cases, and is suppressed as |k| 
is varied. However, in the N — 1 case, there are additional regions of unstable 
q which are not stabilized at these parameter values. 

Figure [7] shows the real part of the Floquet exponents for a second parameter 
set, for which the wave can be stabilized using only one spatially shifted term. 
In this second example, all parameters are the same except that the spatial 
shift Ax! is parallel to the vector (1,0.2), and p\ = 0.1. The two plots show 
the real parts of the exponents of the linearized system with and without the 
temporal feedback (that is, 7 = and 7 = —1/At). There are no positive 
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Qx 



(a) 7 = (b)7 = -l/At 

Fig. 7. The two plots show the real part of the Floquet exponents of the linearized 
system for |k| = 
7 = 0). In [7(E)] 7 



0.284, and N = 1. In 7(a) there is no temporal feedback (i.e. 
= —1/At. The regions with positive growth rate are shaded in, 
as in Figure [U so we can see that with the temporal feedback, the plane wave is 
stabilized. Parameters used are b\ = 2.5, 63 = 2, p% = 0.1, and Axi is parallel to 
(1,0.2). 



growth rates when 7 = —1/At and hence the traveling wave can be stabilized. 



7 Discussion 



In this paper we have examined the possibility of stabilizing traveling waves 
of the CGLE in two spatial dimensions using a combination of both tempo- 
ral and spatial feedback. The feedback is noninvasive, in that it vanishes at 
the targeted wave solution, and is a generalization of that first proposed by 
Pyragas [3] . Our analysis is a local linear stability analysis which involves the 
study of a linear delay equation. The analysis is similar to that in [21], for the 
ID CGLE. However, in the 2D case, the spatial shifts Axj must be carefully 
chosen as otherwise there may be some perturbations q which are unaffected 
by the spatial feedback, meaning the wave cannot be stabilized for any values 
of the control parameters. We show that if the Ax.,- are correctly chosen, the 
results of [21] that the 'best' value of the gain is 7 = —1/At still holds. 

We expect that this method of using noninvasive feedback to stabilize other- 
wise unstable waves could also be applied to amplitude equations other than 
the CGLE. One property of the CGLE which simplifies the analysis is that 
the form of the wave and its dispersion relation can be written down in closed 
form, and are particularly simple. For other amplitude equations this will not 
be the case. A lack of knowledge of the dispersion relation means that it is not 
possible to choose a priori the time delay and two spatial shifts in a consistent 
way so that all temporal and spatial feedback terms vanish simultaneously. It 
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may be possible, however, to circumvent this difficulty if just one spatial term 
and one temporal delay term are employed in the feedback; a single wave may 
then satisfy the criterion that the feedback vanish when it is realized. However, 
the unknown dispersion relation will then be selecting its direction of travel. 

In our analysis, we assumed the feedback parameters 7 and pj were real. This 
is a mathematically convenient choice since the control matrix then ends up 
being a multiple of the identity. However, the choice of real gain coefficients 
for the control terms does not seem natural when all other coefficients in the 
complex Ginzburg-Landau equation are complex. Future work on this problem 
will consider the case of complex gain parameters. 
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